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Abstract. We show that Montecarlo simulations of the TP-AGB stellar population in the LMC 
CN and SMC galaxies using the CB* models produce LF and color distributions that are in closer 
' agreement with observations than those obtained with the BC03 and CB07 models. This is a 
: progress report of work that will be published elsewhere. 
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2 !• Introduction 

C/3 ■ It has been known for quite some time now that intermediate mass stars in the ther- 
mally pulsing asymptotic giant branch (TP-AGB) phase of their evolution contribute 
at least 50% of the NIR light in a simple stellar population (SSP) of age 1-2 Gyr, 
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P\J \ e.g., Maraston (2005) Bruzual (2007, 2011). The treatment of this stellar phase in 
K>" ■ stellar population synthesis models determines the predicted spectral energy distribu- 



tion (SED) of stellar populations in this wavelength and age range. In Fig.[T] we com- 



pare the predictions of Maraston (2005) for a Salpeter IMF, Z = Zq, SSP model at 



ly-j ■ two different ages, with the predictions of three different versions of our code: (a) the 
BG03 models; (b) the minor revision of these models introduced by CB07, and (c) a 
I ^ major revision of this code and models (in preparation, hereafter CB*). The CB07 
■ models use the same sets of stellar tracks and spectral libraries as BC03, except for 
I ' the TP-AGB stars, for which CB07 follow the semi-empirical evolutionary prescrip- 
^ tions by [Marigo fc Glrardi (20()7)] and [Marigo et al. (2008) [ The CB* models used in 



• 'H this paper are based on the stellar evolution models computed by Bertelli et al. (2008), 
/S ■ Tracks are available for metallicities Z = 0.0001, 0.0004, 0.001, 0.002, 0.004, 0.008, 
^ ; Z0=O.O17, 0.04, and 0.07. In CB* the evolution of TP-AGB stars follows a recent pre- 
" " ■ scription by Marigo & Girardi (private communication), which has been calibrated us- 
ing observations of AGB stars in the Magellanic Clouds and nearby galaxies (Girardi 
et al. 2010; [Melbourne et al. 2012| . In the optical range, the CB* models are avail- 
able for the IndoUS (jValdes et al. 2004^ . Miles dSanchez-Blazquez et al. 2006[ ), Stelib 



dLe Borgne et al. 2003[ ), and BaSeL 3.1 (|Westera et al. 2002p spectr al libraries. The NIR 



spectra of TP-AGB stars in CB* are selected from the compilation by Langon & Mouhcine (2002) 



the NASA Infrared Telescope Facility (IRTF) library ( [Rayner et al. 2009[ ), and the C-star 
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model atlas by Aringer et al. (2009) The effects of mass loss and reddening in the spectra 



of TP-AGB stars have been included in these models as described by Gonzalez-Lopezlira et al. (2010) 



The treatment of the TP-AGB in the M05 models is based on the Fuel Consumption 
Theorem and is thus completely independent of the prescriptions used in the BC/CB 
models. 



2. Modeling the LF of TP-AGB stars in the LMC and SMC 

We model the distribution of TP-AGB stars in the CMD in various optical and NIR 
bands for a stellar population of Z = 0.008, close to the LMC metallicity, by means 
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Figure 1. SED in the optical and NIR range for a Salpeter IMF, Z = Zq, SSP at the age 
indicated inside each frame. The different lines clearly visible in the 1 Gyr frame in the 1-2/im 
range represent from bottom to top the models by BC03, CB* (see text), CB07, and M05. 
Whereas the fractional r ange in flux spanned by these models may reach 100% at 1 Gyr, it is 



less than 15% at 12 Gyr. |Kriek et al. (2010)^ Melbourne et al. (2012), and |Zibetti et al. (2012) 
have shown evidence that the treatment ot the TP-AGB stars in the CB07 and M05 models 
overestimates the contribution by TP-AGB stars in the NIR, favoring BC03 and CB*. 
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of Montecarlo simulations (cf. IBruzual 20021 2010). At each time step the mass formed 
in stars is derived from the LMC star formation history ( [Harris fc Zaritsky 2009 1 . The 
stars are distributed in the CMD according to the isochrones computed with the CB*, 
CB07, and BC03 models described in §1. Fig. [2] shows a comparison of the LF derived 
from our three simulations and the observed SAGE data set (jSrinivasan et al. 2009P 
in the IRAC [4.5]/im band. The corresponding IRAC [3.6] — [4.5]/im color distribu- 
tions are shown in Fig.|3l Using the same procedure and the SFH of the SMC from 



Harris & Zaritsky (2004) we model the TP-AGB stellar population in the SMC galaxy, 
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Figure 2. LF derived from our Montecarlo simulations of the star formation history of 
the LMC using the C B*, CB07, and BC03 models, and from the observed SAGE data set 
l|Srinivasan et al. 2009P in the IRAC [4.5]/^m band. The left column corresponds to all the 
TP-AGB stars. In the central and right column only the 0-rich and C-rich TP-AGB stars are 
shown, respectively. For the simulations we assumed the Salpeter IMF, Z = 0.008 isochrones, 
and we kept all the stars with apparent K ^ 12 mag. The heavy gray-line corresponds to the 
SAGE LF. The hea vy black-l ine corr esponds to the simulation LF using the central value of the 
LMC SFH ( [Harris fc Zaritsky 2009[ ). The bracketing light black-lines correspond to the upper 
and lower limit of the SFH derived from the error bars given by these authors. 
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Fig-El In the case of the SMC the chemical evolution indicated by Harris & Zaritsky (2004) 
is included in our simulations. Inspection of Figs. 2-4 show that the LF's computed 
with the CB* models are in closer agreement with the observations that those com- 
puted with the BC03 and CB07 models. These results are consistent with the findings 
by |Kriek et al. (2010)[ Melbourne et al. (2012), and |Zibetti et al. (2012)] (see caption to 
Fig.[T]), and support our choice for the treatment of TP-AGB stars in the CB* models. 
We do not have at hand enough information (isochrones) to perform the same kind of 
comparison with the M05 models. Details of this work will be published in a coming 
paper. 

G. Bruzual acknowledges support from the National Autonomous University of Mexico 
through grants IA102311 and IB102212. 
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Figure 4. Same as Fig.[2]but includ ing the results of the simu lations for the SMC using the 

SFH of|Harris fc Zaritsky (2004)] 
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